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We examine the motion of a photon in the gravitational field of a binary system. The equations 
of motion are geodesic equations in a Schwarzchild background with a tidal force. We specialize the 
equations to that of an edge-on binary and use the method of osculating elements to integrate them. 
This work helps to identify a binary system through the gravitational light deflection of one member 
in the gravitational field of the other member. It is found that the effects of the companion body 
on a photon passing the edge of the star can be potentially detected by astrometric satellites with 
pas precision, if the ratio of the Schwarzchid radius to the star radius, > 10“®. Two different 
cumulative effects on the photon path are identified. 


I. INTRODUCTION 

How is the orbit of a photon that is affected by the 
field of a binary system? This is a three- body problem 
with gravitational interactions. Adding one body to the 
famous two- body problem of a photon in a gravitational 
field brings a remarkable complexity to the problem. It 
is however believed that in a complete treatment of such 
motions and with todays highly accurate astronomical 
observations in hand, these effects must be taken into ac¬ 
count. This subject has been a field of interest since the 
astrometric space mission like GAIA planned to position 
the stars and celestial objects with one microarcsecond 
level of accuracy [T]. At this level, many subtle effects 
are potentially observable. One of them would probably 
be the deflection of light by the field of a binary system. 
When a photon moves in the (weak gravity) field of an 
isolated spherical body, the leading order corrections to 
the photon (Newtonian path) are scaled by where R 
is the closest distance to the gravitating body the photon 
can reach and m is the mass of such body. Considering 
the scales characterizing our problem, it is tempting to 
ask whether there are any relativistic effects produced 
by the second body in the binary that we should be con¬ 
cerned of, or the traces left by the companying body still 
remains a challenge for future missions? 

In this work, we assume that photons move on the null 
paths and are influenced by a binary system, with two 
equal masses, in a circular motion. We also assume that 
the radius of the binary orbit is large, compared to the 
distance between photon and the member of the binary 
which we call it the main body. In adequately weak grav¬ 
ity regimes, effects like the deviation of the photon tra¬ 
jectory from straight line, may be well described by the 
post-Newtonian approximation. The latter assumption 
permits an approximate solution to our problem by re¬ 
ducing the three-body problem to the motion of a photon 
in an effective field of the orbiting bodies. 

To better discuss the dynamics of this system, we intro¬ 
duce two lengthscales. One scale is given by M := 


and the other scale is the radius of the circular orbit, de¬ 
noted by Tq. Let p be the separation between the photon 
and the black hole. We assume that i) M/p is small and 
ii) the second body in the binary is remote. In a refer¬ 
ence frame that moves with the main body, the effects of 
the second body on the geometry of the spacetime can be 
described by tidal potentials. These are given by a lin¬ 
ear perturbation about the Schwarzchild solution; so that 
the perturbed metric satisfies the vacuum field equations, 
perturbatively. 

It is straightforward to apply tidal effects to a system 
comprising a photon and a gravitating body, but solving 
the equations corresponding to the path of the photon 
in this perturbed spacetime would be tricky. A suitable 
framework for solving the path based on osculating or¬ 
bits has already been suggested in [2] . In this framework, 
the motion of the photon is at all times described by 
a family of null geodesics, called osculating trajectories, 
characterized by a set of constants of motion. A nice re¬ 
formulation of the various aspects of the original problem 
can be achieved by the evolution of these constants, if we 
know the interacting forces. The evolution of this set of 
constants is constrained by a set of first order differential 
equations, known as the osculating conditions. In this 
paper we calculate the tidal forces acting on the photon 
employing a metric calculated in a paper by Taylor and 
Poisson [3j. In their work, the main gravitating body 
was a black hole, they used this metric to find the dy¬ 
namics of a black hole, i.e, those produced by an influx 
of gravitational energy across the horizon (like tidal heat¬ 
ing or increase in mass, angular momentum and surface 
area 0) of the black hole. Our interest in this paper is 
mostly in the consequences of the tidal interactions on 
a photon orbital dynamics. The metric given by Taylor 
and Poisson is not necessarily restricted to the weak field 
approximation and can be equally used in small hole ap¬ 
proximation, in which the gravitating body is assumed to 
be much smaller than the second body. In our problem, 
we need to know the tidal fields and their effects on the 
structure of spacetime around the gravitating body, all 
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within the weak field approximation to general relativity. 

Assuming that the photon path would not be too dif¬ 
ferent from the Schwarzchild light paths, we will include 
the force components in the osculating equations given 
in [5] to obtain the observable variations at any point of 
the path. Our results can be evaluated in any coordinate 
systems; this guarantees that those are observable quan¬ 
tities. For example, we find that the total deflection an¬ 
gle depends on the second body position and varies from 

, 3 " 


a minimum equal to 4^ 


1 — TT 




to a maximum 


4^1 

R, 


l + 27r ^ 




where i?* is the closest distance to 


the main body that the photon could reach, if it were an 
isolated body. 

The structure of this paper is as follows: In Section 
II, we start with the metric describing the tidal envi¬ 
ronment around a gravitating body. We find the tidally 
perturbed metric in Schwarzchid coordinates by a coor¬ 
dinate transformation (from harmonic coordinates). We 


would like to focus on the leading effects and neglect the 
smaller ones; so we find the simultaneously weak field 
and slow motion limit of the metric. These perturba¬ 
tions to the metric are considered as static fields on a 
flat background. We then calculate different components 
of the perturbing force and the equations governing the 
photon path in section III. In section IV, the evolution 
of constants of motion is discussed. The discussion will 
be concluded with the summary and conclusions of the 
results in section V. 


II. TIDAL ENVIRONMENT AROUND A 
GRAVITATING BODY 


The general description of the spacetime around a 
spherical gravitating body in an arbitrary tidal environ¬ 
ment is given by a metric, which has the following form 
in Schwarzchild coordinates, (t, p, 9^ = {9, </>)), 
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( 2 . 1 a) 

( 2 . 1 b) 

( 2 . 1 c) 


( 2 . 1 d) 


( 2 . 1 e) 


( 2 . 1 f) 


Here M is gravitational radius of our body, = 

which satisfies = 7 “^ and = 

0. Furthermore, := ^ = 

(sin 0 cos(;li, sin 0 sin(?!),cos 0 ) is a unit radial vector, £Iab 
is the metric on the unit two-sphere, a:“’s are coordinates 
of a Cartesian system, and jab = Sab — £ia^b is a pro¬ 
jection tensor on the direction orthogonal to The 
domain of validity of the metric is the spherical body 
neighborhood. To find this metric, we used the results 
of [3], and performed a coordinate transformation from 
the harmonics coordinates (a:°, a:“) to (t, p, 9, (j)), given by 
= ct,x°’ = { 9 ^^, r = p — M. In [3], Taylor and 

Poisson, examined the motion and the tidal dynamics 


of a black hole placed within a post-Newtonian external 
spacetime. Appendix A presents some crucial results of 
that paper. Our work is a continuation of their effort, 
and this coordinate transformation brings the metric of 
the distorted black hole (A.ll to the new form given in 
(2.1). This metric applies to the spacetime outside any 
spherical body immersed in a tidal environment. It is 
linearized about the Schwarzchild solution and satisfies 
(approximate) vacuum field equations. The tidal envi¬ 
ronment is presented as a functional of arbitrary exter¬ 
nal potentials {£^, £^, B^} built from the lowest order 

tidal moments, £ab and Bab- Within this approximation, 
the metric is accurate to all orders in ^ and the neglected 
terms involve third order in that come from higher or- 
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der tidal moments or the time derivative of quadrapole 
moments. Although the potentials that characterize the 
tidal perturbations are completely general, in this work 
we require that the gravitating body be a member of a 
binary system. 

In studying the path of a photon through a curved ge¬ 
ometry, it is a well known practice to consider the met¬ 
ric perturbations as fields defined on a flat background 
spacetime. In this way the tidal environment is described 
by the weak field limit of the metric. Following the same 
practice, we take ^ as well as the potential terms pro¬ 
duced by the second body, i.e the linear perturbation 
terms to flat Minkowski space, to be small. For the prob¬ 
lem under study, the radius of the binary is large and the 
terms beyond the quadrapole tidal moments were already 
assumed to be small^. For this procedure to make sense, 

we should invoke the slow motion limit, - = \ 1 as 

well. We introduce two dimensionless parameters which 
are almost equally small and denote them by e 



In this limit, the nonzero components of the metric are 
thus written as 

gtt = — p^£'^ + O (e^) , (2.3a) 

= 1 + 2^ - + o , (2.3b) 

gAB = [I - + o (e^)] , (2.3c) 


where and Sab is a symmetric-trace free 

tensor, which its components for circular orbit motion are 
calculated by matching the local metric with the global 
metric that include the gravitating body and the external 
spacetime [3] . These components are given in Appendix. 
By substituting (A.3) in (A.2|, the scalar potential, S'^, 
is obtained as 


£9 = A [l — 3 cos^ 0-1-3 sin^ 9 cos 2 {(j) — iIj) + O (e)] , 

(2.4) 

where A := — ^ is a dimensionful parameter and ip = 
where uj denotes the angular frequency of the tidal 
moments with respect to the frame moving with the main 
body. Before exploring the behavior of the photon in this 
spacetime, there are some points that we would like to 
add: 


^ The light defection in the post-linear gravitational field of two 
bounded masses for the case when the impact parameter is much 
larger (5 times or more) than the distance between binary mem¬ 
bers, is calculated in [5]. In the following, we show that the linear 
terms are so small that we can confidently neglect the post-linear 
terms. 


1. Although the gravitomagnetic tidal (vector) poten¬ 
tials exist in a relativistic description of the geom¬ 
etry around a gravitating body, they do not ap¬ 
pear in weak field limit and slow motion context 
(ri/c <C I, where v is the tidal environment velocity 
scale). This is also true for gravito-electric vector 
and tensor potentials. 

2. The spherical symmetry of the Schwarzchild geom¬ 
etry is perturbationally broken; so the Birkhoff the¬ 
orem will let us to have a non-static vacuum space- 
time. 


3. The tidal geometry is described by a scalar poten¬ 
tial and each scalar potential transforms as such 
under parity. This ensures that the deformed met¬ 
ric be conserved under a parity transformation. 


4. If we expand the tidal potential given in (2.4), in 
terms of spherical harmonics we will have 





-Y° + {Yi + Y^'^) + O (e) 


(2.5) 


It is observed that the coefficients of Y,^^ are zero, 
and the metric components are invariant under 
9 -A TT — 9. This fact demonstrates an additional 
symmetry with respect to the plane 9 = ^. The 
presence of the second body has greatly reduced 
the isometry (Lie) group of Schwarzchild spacetime 
to (discrete) subgroups, namely parity and mirror 
symmetry. 


III. NULL TRAJECTORIES IN THE 
PERTURBED SPACETIME 


The first step before writing the null geodesic equations 
is to calculate the Christoffel symbols for this metric. The 
perturbed geodesic equations will then be written as 


d^xf^ -a dx''dx‘^ 

-pM- 

dA2 dX dX 


(3.1) 


where 


dX dX '■ 


pM = pM I JpM 
2A(T L'O* ' ^ l/a 


(3.2) 


Here Fs and AFs are components of the Schwarzchild 
Christoffel symbols and their perturbations, respectively. 
Appendix B gives the full list of these symbols in the 
weak held limit. The right hand side of (3.1) represents 
the perturbing force components. 

It does not seem to be much hope for solving the set 
of equations (3.1) in a general case. Fortunately, there 


are some special cases in which the remaining symmetries 
simplify the task. It is easy to see that the plane 0 = f 
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satisfies the equation for 0 (A) ? Therefore, if we special¬ 
ize to the case in which photons move in equatorial plane 
of the binary system, i.e, if the inclination parameter, *, 
is exactly 7r/2 and the binary orbit is edge-on, the di¬ 
rection of angular momentum will be conserved and the 
photon will not leave the plane, because F® = 0. Other 
components of the tidal force applied to the photon will 
then be 


F* = + P' - i') 


— 2pip£'^ — + O (e^) 

(3.3a) 

FP = p (jp- -P- gq 


+ p" + pd^.S'^) + O (e^) 

(3.3b) 

F'^ = l (fPpP -fP- P) « 


+ p^dtE^ + 2ppP£‘i + O (e^) . 

(3.3c) 


Even with above simplifications, integrating the equa¬ 
tions of motion and determining the photon path will 
not be simple. Although the perturbing forces are com¬ 
plicated, we believe that these forces keep the form of 
the photon path. Therefore, the leading correction to 
the Schwarzchild null curves can be determined by the 
method of Variation of Parameters (VoP). As described 
in [2], in this method one employs a known solution of 
a fiducial system of differential equations as an ansatz 
for solving the system under study. The known solution, 
Xq, plays the role of instantaneously tangential (or oscu¬ 
lating) trajectories to the true path, (A). Each point 
on the true path with parameter A is coincides geometri¬ 
cally and dynamically with an osculating trajectory with 
one specified set of constants, the number of these 
constants is equal to the number of phase space inde¬ 
pendent variables. The constants are furnished with A 
dependence and evolve along a flow line on a cross sec¬ 
tion of phase space that is consistent with null condition. 
Each and every flow line is constrained by the osculating 
conditions 


ra ^ pM 




dP 


= 0 , 


(3.4) 


but different flow lines may correspond to very different 
interactions. This reformulation of the problem provides 
an insight into the effects of perturbing forces on various 
aspects of the motion. In the framework of VoP, the 
perturbing forces are not assumed to be small and the 
only restriction is that the forces must preserve the shape 
of the path. 

Following [2], these constants are called osculating 
elements, and the Schwarzchild null geodesics family. 


If we calculate in} for = 6, different term s inc lude either 
or which can be easily verified from i B.l I to be zero 
st 6 = 7^12. 


with parameter A, is chosen as Xq. We recall that in 
weak field approximation, the relations describing forced 
Schwarzchild null geodesics are 


( 0 ) = 


sin ((/>-A ((/))) 

R{^) 

3M 

2i?2 (</)) 


1 -f -cos2((/)- A(0)) 


<^(A) = A(A) + 
t(A) = r(A) + 

7^2 

•2 _ p 2 ^ 

P — tj -- 
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(3.5a) 

(3.5b) 


(3.5c) 


(3.5d) 


Here dot denotes derivative with respect to A and 
the osculating elements on these trajectories are 
{F, L, R, A, T}. Some osculating elements, like R, T and 
A, that characterize the initial position on the hducial 
trajectory, are usually referred to as the positional ele¬ 
ments; these elements typically change whenever the ge¬ 
ometry of spacetime is affected with a new interaction. 
Other elements, like L and F, which are called princi¬ 
pal elements, determine on which trajectory the photon 
is moving. We can ignore T and evolve t explicitly by 
using 


t = 


E 

F1- 2—’ 

p 


(3.6) 


which is equivalent to the evolving oit — T, directly. The 
flow lines that satisfy (3.4) are given by 

i ?2 


F = 


A = - 

L = 

R = 


1 

L sin^ {4> — A) 

1 


F2 sin (</) — A) 
i?4 1 


F^ + 0{e^ 


FP + O (e^) 


L sin'* {(p — A) 

i?* COS {(j) — A) 

sin^ {(j) — A) 
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F^ + O (e^) , 
FP + O (e^) , 


(3.7a) 

(3.7b) 

(3.7c) 

(3.7d) 


t = sin {(j) - A) [sin {p - l\) + AM] + O {(?) . 

(3.7e) 


In these equations, ip has been used instead of A (or t), 
as an independent variable and prime denotes derivative 
with respect to p. The equations governing the evolution 
of E and F could have been calculated by the covariant 
formulation of osculating conditions [B]. 

The variations in the osculating elements are assumed 
to be small, therefore one can achieve a good approxima¬ 
tion of these elements by substituting the leading order 
expressions, 

F* 

^ sin((() —A*) 


(l + 0(e)) 


(3.8a) 
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i = E^{l + 0{e))., 




= ^Bi.2(,f-A.) + 0(e) 


(3.8b) 

(3.8c) 

(3.8d) 


into the right hand side of equations (3.71. Here the set of 


constants that the photon trajectory could have assumed, 
if the main body were in a complete isolation, are denoted 
by {i?*, L», i?*. A*}. The different components of the 
perturbing force in the plane 9 = n j2 has the form of 


F* = 

PP = -2E,LU 


i cot (()) — A*) [1 + 3 cos 2 {(j) — tp)] + 3 sin 2 {(j) — Ip) +0 (e^) , 

rC 

k cos {(p — A*) [1 + 3 cos 2{(p — tp)]+3 cos {(p — A*) cos 2{(p — ip) 


P'^ = 2E^A 


—k [1 + 3 cos 2 ((p — Ip)] sin2 (<p — A,) + 3sin2 ((() — ip) 


E^Rluj 


.sim {(p - A*) 


1 + sin^ {p - A*) 


(3.9a) 

— F* + O (e^) , 
(3.9b) 
+ 0(e2). (3.9c) 


Here k := EifR^,/L^, is a dimensionless parameter, which 
its value can be found by substituting, p = i?* in (3.5dl, 
k = 1 + 0 (e). In the derivation of the force components. 


we have used the square root of (3.8c) with a minus sign, 
p = —E* cos {p — A*) + O (e). 

We recall that in preceding equations, p is a function 
of t; so formally, its explicit dependence on p must be 
considered. However, these equations will turn out to be 
simpler by recognizing that in the problem under study, 
osculating elements undergo two types of changes. The 
first is an oscillation with a period equal to the period of 
tidal moments (or multiple of it); the second is a steady 
drift in p. In other words, we have two timescales; t 
and uj~^. On short timescale, t Ki M, the photon moves 
on a geodesy of the Schwarzchild spacetime with a static 
tidal field, characterized by {F, L, i?. A}. Over longer 

timescale, uj~^ « the displacement of the 

V o e S 

second body changes p and causes the path to evolve. 
It is clearly desirable to compute an accurate trajectory 
for the whole domain of the validity of p (or t). How¬ 
ever, the linear perturbation method is limited to pro¬ 
ducing a model for the trajectory over rapid timescale, 
(or an snapshot of the trajectory,) in which the trajectory 
is described by the deviations from the pivot elements, 
{F*, F*, i?*. A*}, and the new configurations (of the sec¬ 
ond body) are neglected. Such approximations fall out 
of phase with the true path after the longer timescale 
« M. In fact a rigorous prescription for the path 
can be provided by a two-timescale analysis, but con¬ 
sidering the order of accuracy by which the short term 
secular changes are captured, at the cost of discarding 
long term effects, it is sensible to assume that p p or 


U)t 


RlE, 


1 


sin {p- A^) 


< 1 . 


(3.10) 


After invoking this approximation, the components of 


the perturbing force can be written as 
F* = 2F*F*Acot {p — A*) 


g Cos {p-2p + A^) ^ 


cos {p — At) 


(3.11a) 


PP = 2F*F*Asin {p — A*) 
F* = 2El A cos {p- At) 


'^ sin{p -2p + At) _ ^ 


sin {p — At) 


(3.11b) 


^ sin {p - 2p + At) _ ^ 


sin (p — At) 


(3.11c) 


The evolution of the conserved quantities will then be 
reproduced by substituting the preceding equations into 
pTl) as 


A = -2ARi 


^ sin {p-2p + At) _ ^ 
sin {p — At) 


(3.12a) 


^ = 2ARl cot {p- At) 


Ft sin^((/) —A*) 


= 21 n2 cos (</>- A^) 
Lt * sin^ ((/) — A*) 


^ sin {p-2p + At) _ ^ 


sin {p — At) 


(3.12b) 


'^ cos {p-2p + At) ^ 
cos {p- At) 

(3.12c) 

P sin{p -2p + At) _ ^ 
sin {p — At) 

(3.12d) 


It is easily seen that the corrections imposed to the 
Schwarzchild null paths by the second body are of or¬ 
der ARl, that is These corrections are so 

small that, even for photons with large impact param¬ 
eters, we can neglect higher order terms in (2.3) with 


confidence. These equations appear to be singular at 
(/) = A*, TT + A*. The singularities correspond to the fact 
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that Schwarzchild geodesics loose their geometric mean¬ 
ing for asymptotic values of tp. 

In the following section, we will give expressions for 
osculating elements; once these functions are known, the 
path is obtained from equations (3.5). Considering the 


long separation of the two timescales, we assumed that 
the geometry around the main body has time transla¬ 
tion symmetry. Before concluding this section, we are 
interested in using this symmetry to calculate energy of 
the perturbed system, E, as a constant on the path. 
The Killing vector whose conserved quantity is energy, 
is K = The conserved energy is then given by 


E = = i[l-‘^+ + O (e^ 


(3.13) 


This is a perturbationally constant quantity even for the 
paths outside of the equatorial plane. Moreover, for the 
trajectories constrained to the equatorial plane, the di¬ 
rection of angular momentum is also conserved. 


IV. EVOLUTION OF OSCULATING 
CONSTANTS 


The expression for the evolution of A in (3.12a I can 
now be used to obtain the explicit dependence on p) at 
each 4>. This gives 


A {P) = A^ + 2ARl cos 2 (V’ - A*)) 

-I- 3 sin 2 (^/i — A*) In (sin |A* — 0|)]. (4.1) 


Fig.[T] shows the variations of A at different points of 
the path for different ip directions at the initial condition 
A* = 0. The contribution of the second body in the 
total deflection angle of a light ray passing near a binary 
system is 


^ A*+7r 


Ad(j) 


2ttARI [1 — 3 cos 2{'tp — A*)]. 


(4.2) 


This implies that for edge-on binaries, the second body 
gives maximum increase in the deflection angle when 
Ip — A^ = {f,^} and maximum decrease is obtained 
when i/) — A* = {0, tt}. In [7], it is shown that among or¬ 
bits with arbitrary inclination, the maximal possible light 
deflection is attained for edge-on binaries, (i = tt/2). Our 
calculations complement the criteria for orbital elements 
by using those the maximal changes in positional ele¬ 
ments are determined. The various gravitational effects 
in the light propagation are additive in weak field limit; 
so for the total deflection angle, <5, we will find that 


4^ - AnARl <S< 4^ -f SttARI (4.3) 


By integrating (3.12b), one can obtain 



FIG. 1: Variations of A —A, vs. (p and p}. Here A* is assumed 
to be zero. On the vertical axis, we have chosen ARI = 1 



FIG. 2: Variations of R/R* — 1 vs. (p and ip. On the vertical 
axis, we have chosen ARii = 1 


— (1 — 3 COS 2 (A* — tp)) Insin^ {(p — A*)] . (4.4) 
This amounts to a total change in R equal to 

127rAi?3sin2(A* - V'). (4.5) 


The photon gets a maximum value of total R- deflection 
when Ip — A^ = change in total R is 

expected when '(/) — A* = {|,7 r,^}. 

Obviously, one can find a different parametrization of 
the path by first integrating ( |3.7e| ) and keeping the lead¬ 
ing order terms to produce. 


t{(P) 


Ej^ 

L,Rl 


1 

4 


sin 2{(p — A*) -I- 4M cos {(p 


A*) 

(4.6) 


R 

R* 


1 -|- 2ARI [6 sin2 (A* — ?/>) (1 — cot {cp — A*)) 
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FIG. 3: Variations of — 1 vs. (p and tp. On the vertical FIG. 4: Variations of L/L* — 1 vs. cp and -ip. On the vertical 

axis, we have chosen AR^ = 1 axis, we have chosen AR^ = 1 


and then omitting (p in favor of t, numerically. The ex¬ 
pressions describing the evolution of other osculating el¬ 
ements, up to the leading order in e, are 


E 

E^ 


L 

L:^ 


1 -h ARI [6sin2 (A* — ip) cot {(p — A*) 


1 -I- 3 cos 2 (A* — Ip) 
sin^ (0 — A*) 


(4.7a) 


1 -I- ARI [—6 sin 2 (A* — ip) cot^ (^ — A,) 
1 — 3 cos 2 (A* — Ip) 


sin^ (0 - A*) 


(4.7b) 


In Figs.[^|^ we show the evolution of the path under the 
influence of the second body and initial condition A* = 0. 
To avoid the singularities at asymptotes, (p = 0, tt, the (p 
variations does not cover the whole domain, also we have 
chosen AR^ = 1. The four figures show the two posi¬ 
tional elements, (A,i?), and the two principal elements, 
{E,L), as functions of (p and ip. The comparison of be¬ 
havior of these elements at asymptotes shows that there 
is a secular change in (A, R); but there is no changes in 
the values of the principal elements in these limits. For 
particles moving on timelike bound orbits, it is believed 
that the changes in principal elements occur when per¬ 
turbing forces include a dissipative part [8j. Likewise, 
equations (4.7b) and (4.7a) imply that L (0) = L (tt) and 
E {0) = E (tt) for massless particles; even though the to¬ 
tal angular momentum is not conserved in our problem. 
Although this result is obtained for a tidal interaction 
and edge-on binaries, it seems to be also valid for general 
inclinations. 


V. SUMMARY AND CONCLUSIONS 


In this paper, we considered a photon in the gravi¬ 
tational field of a binary system and found the form of 
trajectories invoking the weak field and slow motion ap¬ 
proximations. In practice, we found the magnitude of 
secular changes to the path of a light ray propagating in 
the field of a spherical gravitating body, that arises in 
the presence of a companion body. This work illustrates 
the application of the osculating elements formalism to 
the computation of photons paths evolutions in the the 
Schwarzchild spacetime, as well as the tidal geometry 
around a black hole in reducing a photon-binary problem 
to an effective two body problem. This work shows how 
a companion star can affect the curvature of a light ray 
brought about by the main lens. It is known that edge-on 
binaries deflect the light more than others [7]; this cri¬ 
terion is further developed in this paper. Provided that 
the companion body is in a right distance and at a right 
orbital phase, in principle, it can be traced with a sensi¬ 
tive detector. The right distance. To, and phase, ?/> — A*, 
are determined by the geometry of spacetime around the 
main lens. Equation (3.12) indicates that edge-on bina¬ 
ries with To = O JiWj, and ip — = {f,^} can 

produce deflections of order Therefore, these 


corrections have potentially detectable effects on GAIA 
(or any other fxas « I0“® astrometric mission) observ¬ 
ables, if e > 10“®. This condition holds for known typi¬ 
cal binary systems. In this estimation, it is assumed that 
the photon passes by the edge of the star. 

The importance of determining such corrections to 
Schwarzchild null geodesics can be understood when as¬ 
trometric information is read off the light paths at a high 
precision level. The astrometric observations of a GAIA- 


















like satelite has been modeled within the PPN formal¬ 
ism to Post-Newtonian gravity. The sensitivity of such 
astronomy to some PPN parameters m as well as the 
light bending due to quadrapole moment of an axisym- 
metric planet (like Jupiter) [T^] has been established. In 
such feasibility studies, one must consider the light de¬ 
flection in the vicinity of a star which is i) a member 
of a binary system and ii) crosses one of the astromet¬ 
ric fields during the GAIA’s mission. A comparision be¬ 
tween the details of the background field corresponding 
to the time when the star is at the center of the field 
with that when it is not, will help to extract the contri¬ 
bution of the binary system. The difference between the 
two patterns may then be fitted to the deflection model 
to determine whether the effect is detectable or not. In 
spite of the common belief that 50 percent of stars are 
a member of a binary (or multiple) system |I3j and the 
order of correction terms is in GAIA’s accuracy level, 
finding a favourable field and binary system with known 
ephemerides is a big challenge that is out of the scope of 
this paper. On the other hand, the trace of wide binary 
lenses studied in this paper may be found in microlensing 
light curves. If repeating microlensing peaks m appear 
in the microlensing light curves during the GAIA’s mis¬ 
sion period, it will be a good signal for a binary candidate 
with M Tq. Analytical results derived in this paper 
can be used in models htted to the microlensing light 


curves. Such models provide a convenient framework for 
mutual verification of the results in an astrometric field 
common in different sky missions. 


Appendix A TIDAL ENVIRONMENT AROUND 
THE MAIN BODY 

In this appendix we list some results presented in [5] 
which investigated a tidally deformed spacetime around 
a nonrotating black hole. The tidal environment around 
a slowly rotating black hole has been studied recently 
by Poisson To capture the effects associated with 
the dragging of inertial frames, this calculation required 
to go beyond first order perturbation theory. When the 
rotational deformations are applied to a circular binary, 
the tidal moments will incorporate terms that involve 
powers of v/c. Since these two metrics are the same at 
the orders that we require, the higher order terms do not 
affect on the results of our work. Some notations we use 
in our work, and are adapted from [5] , will be introduced 
here. We will express the results in the form required for 
our specific purpose. 

The metric components of the black hole immersed in 
an arbitrary tidal environment, in harmonic coordinates 
(a;°,a;“), are given by 
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(A.la) 
(A. lb) 

(A.lc) 


Here M is the gravitational radius of the black hole, when 
it is in complete isolation, are a system of Gartesian 
coordinates, := ^ = (sin 0 cos ((), sin 0 sin cos 0) is a 
vector which satisfies = 1 and 'jab = Sab — Haf^b 

is a projection tensor on the direction orthogonal to 1/“. 
The contribution of the tidal environment is expressed 
as an expansion in the strength of the tidal coupling. 
The expansion parameter is the inverse of the radius of 
curvature, of the external spacetime in which the 

black hole moves. In the neglected terms, C is the scale 
of spatial inhomogeneity in the external spacetime. At 
the lowest order in TZ~^, the metric can be described by 
a set of external potentials {£'^,built from 


quadrapole moments, £ab and Bab by 

£■« = (A.2a) 

£1 = (A.2b) 

C=^(27aVfcd+7ab£^), (A.2 c) 

(A.2d) 


The tidal moments are formally defined in terms of the 
Weyl tensor of the external spacetime. In practice, 
however, for every application those are determined by 
matching the local metric to the global metric that in- 









9 


eludes the black hole and the external spacetime. For 
the specialized case that the black hole is part of a bi¬ 
nary system in circular motion and the global spacetime 
is expressed in the post-Newtonian gravity regime, such 
a matching were carried out in [S] and m- For the prob¬ 
lem under study in our paper where the two gravitational 
bodies in the binary system have equal masses and the 
relative orbital velocity, u, within the binary system is 

slow, ^ nonvanishing compo¬ 

nents of the tidal moments are given by 
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(A.3a) 

(A.3b) 

(A.3c) 

(A.3d) 


B 23 = B 31 = - '^v sin V’ -b O (c , 
Bis = B 32 = - cos V’ + O (c“^) . 


(A.3e) 

(A.3f) 


Here V’ is the phase related to the angular velocity, w, 
of the tidal moments with respect to the frame moving 

with the black hole, ip = wt = -b O (c)^) ■ 

can also be defined in terms of the orbital velocity of the 
second body, w, and the precessional angular frequency of 
the black hole frame relative to the barycentric frame, H, 
hy ip = (jji — Ht, in which t is the global time coordinate. 


Appendix B CHRISTOFFEL SYMBOLS 

The nonvanishing Christoffel symbols built from met¬ 
ric (2.3) are 


Bit = -r^p = 

np = f 

pt _ 1 „2 d 
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